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Abstract 

Aims. We want to derive the mass-metallicity relation of star-forming galaxies up to z ~ 0.9, using data from the VIMOS VLT Deep Survey. The 
mass-metallicity relation is commonly understood as the relation between the stellar mass and the gas-phase oxygen abundance. 

Methods. Automatic measurement of emission-line fluxes and equivalent widths have been performed on the full spectroscopic sample 
of the VIMOS VLT Deep Survey. This sample is divided into two sub-samples depending on the apparent magnitude selection: wide (/ab < 22.5) 
and deep (/^b < 24). These two samples span two different ranges of stellar masses. Emission-line galaxies have been separated into star- 
forming galaxies and active galactic nuclei using emission line ratios. For the star-forming galaxies the emission line ratios have also been 
used to estimate gas-phase oxygen abundance, using empirical calibrations renormalized in order to give consistent results at low and high 
redshifts. The stellar masses have been estimated by fitting the whole spectral energy distributions with a set of stellar population synthesis models. 

Results. We assume at first order that the shape of the mass-metallicity relation remains constant with redshift. Then we find a stronger 
metallicity evolution in the wide sample as compared to the deep sample. We thus conclude that the mass-metallicity relation is flatter at higher 
redshift. At z ~ 0.77, galaxies at 10^ "* solar masses have —0.18 dex lower metallicities than galaxies of similar masses in the local universe, 
while galaxies at lO'"-^ solar masses have —0.28 dex lower metallicities. By comparing the mass-metallicity and luminosity-metallicity relations, 
we also find an evolution in mass-to-light ratios: galaxies at higher redshifts being more active. The observed flattening of the mass-metallicity 
relation at high redshift is analyzed as an evidence in favor of the open-closed model. 

Key words, galaxies: evolution - galaxies: fundamental parameters - galaxies: abundances - galaxies: starburst 



1. Introduction 

\ The stellar mass and the gas-phase metallicity of a galaxy are 
two of the main parameters involved in the study of galaxy for- 
mation and evolution. As cosmological time progresses, theory 

, predicts that both the mean metallicity and stellar mass of galax- 
ies increase with age as galaxies undergo chemical enrichment 
and grow through merging processes. At any given epoch, the 
accumulated history of star formation, gas inflows and outflows, 
affects a galaxy mass and its metallicity. Hence one expects the- 
ses quantities to show some correlation and this will provide cru- 
cial information about the physical processes that govern galaxy 
formation. 

Fi rst discovered for irregular galaxies (iLequeux et al.l 
Il979l) . th e mass-metalli cit y rel ati on has been int e nsively 
studied (ISkiUman et al.l 1 19891: iBrodie & Huciiral Il99ll: 

* based on data obtained with the European Southern Observatory 
Very Large Telescope, Paranal, Chile, program 070.A-9007(A), and on 
data obtained at the Canada-France-Hawaii Telescope, operated by the 
CNRS in France, CNRC in Canada and the University of Hawaii. 



Zaritskv et al.l 119941: [Richer & McCalll 119951; iGarnett et al.l 
19971: iPilvudn & Ferrinill2000l among others) a nd is now well 
establi shed in the local univ erse by the work of iTremonti et al.l 
( |2004 with SDSS data and iLamareille et al.l ( |2004|) with 2dF- 
GRS data, the latter done on the luminosity-metallicity relation 
which is easier to derive when small number of photometric 
bands are available. These two studies have shown in two 
different ways that the mass-metallicity relation is mainly driven 
by the decrea s e of m etal loss when stellar mass increases. 
ITremonti et alj ( |2004|) have indeed observed an increase of 
the ef fective yield with stellar mass, while ILamareille et al.l 
(12004 have shown the increase of the slope of the luminosity- 
metallicity relation. This last trend ha s also been obs erved 
down to much lower galaxy masses by iLee et al.l (l2006l) . We 
nevertheless note that they also observe a large scatter in the 
effective yield, which they find difficult to analyse in the context 
of more efficient mass loss among low mass galaxies. 

Hierarchical galaxy formation models, that take into account 
the chemical evolution and feedback processes, are able to repro- 
duce the observed mass-metallicity relation in the local universe 



(e.g.lDe Lucia et al.ll2004tlde Rossi et al.ll2007HF"inlator & Davi 
l2008h . However these models rely on free parameters, such as 
feedback efficiency, which are not yet well constrained by ob- 
servations. Alternative scenarios have been proposed to explain 
the mass-metallicity relation including low star formation ef- 
ficiency in low-mas s galaxies caused by supernova feedback 
(iBrooks et alJ l2007h and a variable stellar initial mass func- 
tion being more top-heavy in galaxies with high er star forma- 
tion r ates, thereby producing higher metal yields (iKoppen et al.l 
I2007h . 

The evolution of the mass-metallicity relation on cosmo- 
logical timescales is now predicted by semi-analytic models of 
galaxy formation, that include chemical hydrod ynamic simula- 
tions within the standard A-CDM framework dPe Lucia et alj 
120041: iDave & Oppenheim"e3l2007l) . Reliable observational es- 
timates of the mass-metallicity relation of galaxies at different 
epochs (and hence different redshifts) may thus provide impor- 
tant constraints on galaxy evolution scenarios. Estimates of the 
mass-metallicity - or luminosity-metallicity - rel ation of galax- 
ies up to z 1 are l i mited so far to small samp les ( Hammer et al 



2005 



2003 



Liang et al 



iMaier et al 



2004 iMaier et 



nai l samp les 
aP 12004 E 



obulnickv et al 



I2005L among others). Recent studies have 



been performed on larger sam ples (> 100 galaxies) but with 
contradictory results. Savaglio e t al.l (|2005|) concluded with a 
steeper slope in the distant universe, interpreting these results 
in the framework of th e closed-box model. On the contrary, 
iLamareille et al.l (l2006al) did not find any significant evolution 
of the slope of the luminosity-metallicity relation, while the av- 
erage metallicity at z « 0.9 is lowered by 0.55 dex at a given 
luminosit y, and by 0.2 8 dex after correction for luminosity 
evolution. IShapley et aH (l2005l) and iLiu et all (l2008h have also 
found 0.2 — 0.3 dex lower metallicitie s at z = lin the DEEP2 
sample. At higher redshifts lErb et al.l ( l2006l) derived a mass- 
metallicity relation at z ~ 2 lowered by 0.3 dex in metallicity 
compared wit h the local estimate, a trend which could extend 
up to z ~ 3.3 dMaiolino et al.ll2007l) . Please note that the above 
numbers are given after the metallicities have been renormal- 
ized to the same referenc e caUbration in order to be comparable 
(iKewlev & Ellisonll2008h . 

In this work, we present the first attempt to derive the mass- 
metallicity relation at different epochs up to z 1 using a unique, 
large (almost 20 000 galaxies) and homogeneous sample of 
galaxies selected in the VIMOS VLT Deep Survey (VVDS). The 
paper is organized as follows: the sample selection is described 
in Sect. |2l the estimation of stellar masses and magnitudes from 
SED fitting are described in Sect. |3] the estimation of metal- 
licities from line ratios is described in Sec. IH and we finally 
study the luminosity-metallicity (Sect. [5]), and mass-metallicity 
(Sect.|6]l relations. Throughout this paper we normalize the de- 
rived stellar masses, and absolute magnitudes, with the stan- 
dard A-CDM cosmo logy, i.e. /? = 0.7, i2ni = 0.3 and D.a = 0.7 
dSpergel et al.ll2003h . 

2. Description of the sample 

The VIMOS VLT Deep Survey (VVDS. lEeFevre et al.ll2003bl) 
is one of the widest and deepest spectrophotometric surveys 
of distant galaxies with a mean redshift of z ~ 0.7. The 
optical spectroscopic data, obtai ned wi th the Visible Multi- 
Object Sectrograph (VIMOS, iLe Fevre et al . 2003a) installed at 
ESO/VLT (UT3), offers a great opportunity to study the evo- 
lution of the mass-metallicity relation on a statistically signif- 
icant sample up to z ~ 0.9. This limitation is imposed, in the 
current study, by the bluest emission lines needed to compute a 
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Figure 1. Redshift distributions of our wide and deep samples of 
star-forming galaxies, and of their union (see text for details). 



metallicity (i.e. [Oil] A 3727, Hj3 and [Olll]A5007) being red- 
shifted out of the wavelength range of the survey (approxi- 
mately 5500A < A < 9500A). Metallicities measurements up to 
z w 1 .24 using a different set of lines will be provided in a subse- 
quent paper (Perez-Montero et al. 2008). The spectra have been 
taken in one arcsecond width slits, under a spectral resolution 
« 230. 

The first epoch VVDS spectroscopic sample is purely ap- 
parent magnitude-selected and is divided into t wo deep fields 
(17.5 < /ab < 24) : VVDS-02h (herefat er F02, iLe Fevre et all 
120051) and CDFS dLe Fevre et al.ll2004 : and three wide fields 
(17.5 < /ab < 22.5): VVDS-22h, VVDS-l Oh and V VDS-14h 
(hereafter F22, FIO, and F14 respectively, iGarilU e t al. 2008]). 
Spectroscopic data hav e been reduced using the VIPGI pipeline 
dScodeggio et alJl2005l) . which performs automatic ID spectra 
extraction, correction for telluric absorption lines, and flux cal- 
ibration. We only kept objects with a redshift known at a confi- 
dence level greater than 75% (i.e. VVDS redshift flags 9, 2, 3 and 
4). Duplicated observations were not used (some objects have 
been observed twice or more, either by chance or intentionally), 
we always kept the main observation. We also removed, for the 
purpose of our study (i.e. star-forming galaxies), all stars (zero 
redshift) and the sample of broad-line A ctive Galactic Nuclei 
(hereafter AGNs; iGavignaud et al.ll2006l) . 

The photometric coverage of the VVDS spectroscopic 
sample is as follows: CFH12k obse rvations in BVRI bands 
for F02, F22 FIO and F14 fields (iMcCracken et all 120031: 
iLe Fevre et al. 2004j). completed b v U band observations in 
the F02 field dRadovich et al]l2004 . CFHTLS (Canada-France- 
Hawai Telescope Legacy Survey) u * g'r'i'z' bands for F02 
and F22 fields, and JKs photometry a vailable in F02 field 
(llovino et al.l 120051: iTemporin et al.ll2008l) For the CDFS field, 
we use CFH12k UBV RI dArnouts et al.l l200ll) and HST bviz 
(iGiavaUsco et al.ll2004 observations. Finally, objects in F02 and 
F22 fields have been cross-matched with the UKIDSS public cat- 



alog (IWarren et al.ll2007h . providing additional observations in 
JK bands. 

Throughout this paper, we will use a deep and a wide sam- 
ples. The deep sample is made of F02 and CDFS fields. Up 
io z < 1.4 (the limit for the [Oll]A3727 emission line to be in 
the observed wavelength range), the deep sample contains 7404 
galaxies (non-broad-line AGN) with a redshift measured at a 
confidence level greater than 75%. The wide sample is made 
of F22, FIO, F14, F02 and CDFS fields (the last two being Hm- 
ited out to /ab < 22.5). Up to z < 1.4 the wide sample contains 
13978 galaxies (non-broad-line AGN) with a redshift measured 
at a confidence level greater than 75%. We draw the reader's 
attention on the fact that the two samples overlap for galaxies 
observed in F02 or CDFS fields at /ab < 22.5, giving a total 
number of 18648 galaxies. Fig.[T]shows the redshift distribution 
of the star-forming galaxies (see Sect. 12.21 ) in the wide and deep 
samples and in their union. 

2.1. Automatic spectral measurements 

The emission lines fluxes and equivalent widths, in all VVDS 
galaxy (non-broad-line AGN) spectra, have been measured 
with the platefit_vimos pipeline. Origin ally developed for the 
high spectral resolution SDSS spectra dTremonti et al.l l2004t 
iBrinchmann et al.l |2004|) . the platefit software has first been 
adapted to fit accurately all emission lines after removing the 
stellar continuum and absorpti on lines from lower resol ution and 
lower signal-to-noise spectra (iLamareille et al.ir2006bl) . Finally, 
other improvements have been made to the new platefit_vimos 
pipeline thanks to tests performed on the VVDS and zCOSMOS 
(iLiUv et al.ll2006l) spectroscopic samples. A full discussion of 
platefit _vimos pipeline will be presented in Lamareille et al. (in 
preparation) but, for the benefit of the reader, we outline some of 
the main features in this section. 

The stellar component of the spectra is fitted as a non- 
negative linear combination of 30 single stellar population tem- 
plates with different ages (0.005, 0.025, 0.10, 0.29, 0.64, 0.90, 
1.4, 2.5, 5 and llGyr) and metallicities (0 2, 1 and 2.5 Zpi). 
These templates have been derived using the iBruzual & Chariot! 
(l2003h library and resampled to the velocity dispersion of VVDS 
spectra. The dust attenuation in the stellar population model is 
left as a free parameter Foreground dust attenuation from the 
Milky Way has been corrected using Schlegel et al. (1998) maps. 

After removal of the stellar component, the emission lines 
are fitted together as a single nebular spectrum made of a sum 
of gaussians at specified wavelengths. All emission lines are 
tied to have the same width, with exception of the [Oll]A3727 
line which is actually a doublet of two lines at 3726 and 
3729 A and appear broadened compared to the other single 
lines. The spectral resolution is also too low to clearly separate 
rN lllA6584 and Ha emissio n lines. It has been shown however 
bv ILamareille et al.l (l2006bl) that the [Nii]A6584/Ha emission- 
line ratio, which is used as a metallicity calibrator, can be reli- 
ably measured above a sufficient signal-to-noise ratio even at the 
resolution of VVDS. 

Note that because of the limited observed wavelength cover- 
age of the spectra, we do not observe all well-known optical lines 
at all redshifts. Thanks to the stellar-part subtraction, no correc- 
tion for underlying absorption has to be applied to the Balmer 
emission lines. 

The error spectrum which is needed for both fits of the stellar 
or nebular components is calculated as follows: a first guess is 
obtained from photons statistics and sky subtraction and is calcu- 



Table 1. Definition of the H5w index 
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Figure 2. Comparison between the standard Lick H5a index and 
the new broad index defined for low resolution VVDS spectra, 
namely H5w (see Tab.[T|), for a set of model templates taken from 
iBruzual & Charloa ( l2003l) library and covering a wide range of 
ages and metallicities. The solid curve is the least square fit, 
which foUows Eq.[T] Values are given in A. 



lated directly by the VIPGI pipeline. A fit of the stellar and neb- 
ular components is performed with platefit_vimos using this first 
error spectrum. The residuals of this first fit are then smoothed 
and added to the error spectrum quadratically, and a new fit is 
performed with platefit_vimos. 

Finally, detected emission lines may also be removed from 
the original spectrum in order to get the observed stellar spec- 
trum and measure indices from it, as well as emission-line equiv- 
alent widths. The underlying continuum is indeed obtained by 
smoothing the stellar spectrum. Then, equivalent widths are 
measured as direct integration over a 5 Gaussian-sigma bandpass 
of the emission-line Gaussian model divided by the underlying 
continuum. 

For the absorption lines, one has to be aware that the res- 
olution of VVDS spectra is too low in order to measure accu- 
rately Lick indices: the bandpasses of these indices are narrower 
than the width of the absorption lines at this resolution. In order 
to measure the absorption component in the H5 line which is 
used to derive stellar masses (see below), we have defined a new 
broadened index for this specific line called H5w- Table[T|gives 
the three bandpasses of this index. Figure |2] shows the compar- 
ison between this index and the standard Li ck H^a index mea- 
sured on various model templates taken from lBruzual & Charloll 
(l2003h library and covering a wide range of ages and metallic- 



Table 2. Sets of emission lines which are required to have S/N > 
4 in various redshift ranges associated to the various diagnostics 
used in this study. 



Redshift 



"[On] H^B [Om] Ha \Nu] [Suf 



0.0<z<0.2 
0.2<z<0.4 
0.4<z<0.5 
0.5 <z<0.9 



/ 
/ 
/ 



/ 
/ 
/ 



/ 
/ 



/ 



/ 



* For the red diagnostic the [Nll]A6584 and [Sll]A16717+6731 emis- 
sion Unes may not be used at the same time. 
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Figure 4. Blue spectral classification of 1060 narrow emission- 
line galaxies in the redshift bin 0.5 < z < 0.9. The emission-line 
ratios - i.e. [Olll]A5007/Hj3 and [Oll]A3727/Hj8 - are calcu- 
lated using equivalent wid ths. The red solid curve is the empir- 
ical separation defined by iLamareille et al] ( |2004|) . The dashed 
curves delimits the error domain where both star-forming and 
AGNs galaxies show similar blue line ratios. The star-forming 
galaxies are plotted as blue triangles, the Seyfert 2 galaxies as 
green solid squares, the candidate star-forming galaxies as red 
solid pentagons, and the candidate AGNs as orange open pen- 
tagons. The error bars are shown in grey. 



ities. For future comparison purpose between VVDS measure- 
ments and other studies, we derive the following relation: 



H5a = -4.69+1.691 xH5w- 0.0447 11 xH5^ 



(1) 



Note that the majority of points lying outside the fitted re- 
lation in Fig. |2] (for H5w < 5A) are models with a very young 
stellar population (< 10 Myr), hence this relation is only valid 
for stellar populations older than 10 Myr. 

2.2. Selection of star-forming galaxies 

2.2.1. Description of the various diagnostics 

To ensure accurate abundance determinations, we restrict our 
attention to galaxies with emission lines detected at S / N > 4, 
where the set of lines considered for this S/N cut varies with 



redshift. Table |2] provides a summary of the emission lines re- 
quired to have S/N > 4 in the various redshift ranges associated 
to the various diagnostics defined below. 

Our work intends to understand the star-formation process in 
galaxies. Thus we have to remove from the sample of emission- 
line galaxies those for which the source of ionized gas, respon- 
sible for these lines, is not hot young stars. Emission-line galax- 
ies can be classified in various spectral types, which depend on 
the nature of their nebular spectrum. The two main categories 
are the star-forming and AGN galaxies. Their source of ionizing 
photons are respectively hot young stars, or the accretion disk 
around a massive black-hole. As stated before, broad-line AGNs 
(also called Seyfert 1 galaxies) are already taken out of our sam- 
ple. The narrow-line AGNs are divided in two sub-categories: 
Seyfert 2 galaxies , and LINERs . 

Both Seyfert 2 galaxies and LINERs can be distinguished 
from star-forming ga l axies using standard diagnos tic diagrams 
(iBaldwinet al.lll98ll: IVeilleux & Osterbrockl [l987l) . which are 
based on emission-line ratios. The most commonly used classi- 
fication is provided by the [Oiii]A5007/Hj3 vs. [Nii]A6584/Ha 
diagram (Fig. [3^), for which a semi-empirical limit between 
star-form ing, Seyfert 2 galax ies, and LINERs has been de- 
rived by iKewley et al.l (1200 ll) from photo-ionization models. 
Due to our low spectral resolution, the [Nll]A6584 line is not 
always detected and its detection is less reliable than other 
lines, because of possible deblending problems with the bright 
Ha line. Consequently, most of the galaxies in the same red- 
shift range are actually classified with the [Olll]A5007/Hj3 vs. 
[Sii]AA6717H-6731/Ha diagram (see Fig.[3b). 

The two diagrams shown in Fig. [3] which we also refer to 
as the red diagnostics, can only be used with our data in the 
0.2 <z< 0.4 redshift range, as the desired emission lines are not 
visible outside this domain. At higher redshifts, one can use the 
alternative blue diagno stic which was empirically calibrated by 
ILamareille et al.l (|2004|) from 2dFGRS data, in order to give sim- 
ilar results as the red diagnostics. The blue diagnostic is based 
on the [Olll]A5007/Hj3 vs. [Oll]A3727/Hj8 diagram (line ratios 
calculated on rest-frame equivalent widths), and may be applied 
in the 0.5 < z < 0.9 redshift range (see Fig.|4]i. 

For the other redshift ranges, where neither the red nor the 
blue diagnostic applies, we use a minimal classification based on 
Ha, [Nii]A6584 and [Sii]AA6717h-6731 for 0.0 < z < 0.2, and 
on Hj3 and [Olll]A5007 for 0.4 < z < 0.5. We propose the fol- 
lowing selection (called H)3 diagnostic) for star-forming galax- 
ies in the 0.4 < z < 0.5 redshift range: log([Oiii]A5007/Hj3) < 
0.6 (rest-frame equivalent widths, see Fig.|5h). We note that real 
star-forming galaxies which are lost with the Hj3 diagnostic are 
mainly low metallicity ones. In any case a quick check on SDSS 
data has told us that no more than 40% of star-forming galax- 
ies with 12 + log(0/H) < 8.1 actually fall in the AGN region 
of the H)3 diagnostic, this proportion being negligible at higher 
metal lici ties. 

Figure [Sji shows the minimum classification for the 0.0 < 
z < 0.2 redshift range (called Ha diagnostic). The pr oposed sep- 
aratio n has been derived based on the 2dFGRS data dLamareilld 
|2006|) and follows the equation (rest-frame equivalent widths): 

log([Nii]A6584/Ha) = 

-0.4 if log([Sii]AA6717- 
-0.7-1.051og([Sii]AA6717 



6731/Ha)> -0.3 
-6731 /Ha) otw. 



(2) 



This equation have been derived to efficiently reduce both 
the contamination by real AGNs in the star-forming galaxies re- 
gion (< 1 %), and the fraction of real star-forming galaxies which 
are lost (< 4%). 
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(a) (b) 

Figure 3. Red spectral classification of 412 narrow emission-line galaxies in the redshift bin 0.2 < z < 0.4. The emission-line ratios 
- i.e. [Oiii]A5007/H)3 and [Nii]A6584/Ha (a) or [Oiii]A5007/Hj3 and [Sii]AA6717+673 1/Ha(b) - are calculated using equivalent 
widths. The red solid curves are the s emi-empirical separations defined by iKewley et al.l (1200 ih . The dashed line is the separation 
proposed by lKauffmann et aTl (12003 al) . The star-forming galaxies are plotted as blue triangles, the Seyfert 2 galaxies as green solid 
squares, the LINERs as green open squares. The error bars are shown in grey. 
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(a) (b) 

Figures. Intermediate spectral classifications of 1479 (a) and 43 (b) narrow emission-line galaxies. The two plots show the inter- 
mediate diagnostic diagrams, used in the redshift bins 0.4 < z < 0.5 (a) and 0.0 < z < 0.2 (b), or in other redshifts bins for galaxies 
with insufficient number of detected emission lines (therefore not classified with the red or blue diagnostic). The emission-line ratios 
- i.e. [Olll]A5007/H)3, fSlllAA6717-H6 731/Ha, and [NllU 6584/Ha - are calculated using equivalent widths. The red solid Hnes 
are the empirical separations defined bv iLamareilld ( |2006|) . In the left panel, the dashed lines delimits the error domain where both 
star-forming and AGNs galaxies show similar blue line ratios. The star-forming galaxies are plotted as blue triangles, the Seyfert 2 
galaxies as green solid squares, the candidate star-forming galaxies as red solid pentagons, and the candidate AGNs as orange open 
pentagons. The error bars are shown in grey. 



Table 3. Statistics of star-forming galaxies and narrow-line AGNs among emission-line galaxies for various diagnostics depending 
on which emission lines are observed (see text for details). The results are presented for the two wide and deep samples used in 
this paper, and for their union (some objects are in common). We also mention for each sample the total number of objects which 
includes emission-line, faint and early-type galaxies. 



Sample/Diagnostic 


red 


% 


blue 


% 


Ha 


% 




% 


all 


% 


wide (13978) 






















emission-line (total) 


455 


100 


924 


100 


47 


100 


1402 


100 


2828 


100 


star-forming 


415 


91 


768 


83 


30 


64 


978 


70 


2191 


77 


candidate s.-f. 








103 


11 








306 


22 


409 


14 


candidate AGN 








35 


4 








74 


5 


109 


4 


Seyfert 2 


18 


4 


18 


2 


17 


36 


44 


3 


97 


3 


LINER 


22 


5 




















22 


1 


deep (7404) 






















emission-line (total) 


151 


100 


568 


100 


11 


100 


574 


100 


1304 


100 


star-forming 


136 


90 


412 


73 


10 


91 


312 


54 


870 


67 


candidate s.-f. 








92 


16 








166 


29 


258 


20 


candidate AGN 








47 


8 








47 


8 


94 


7 


Seyfert 2 


5 


3 


17 


3 


1 


9 


49 


9 


72 


6 


LINER 


10 


7 




















10 


1 


union of wide and deep (18648) 






















emission-line (total) 


469 


100 


1213 


100 
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The results of all diagnostics, in the wide, deep, and global 
samples, are shown in Table. [3] 



that the contamination of star-forming galaxies by LINERs, in 
the blue or Hj3 diagnostic, is less than 1%. 



2.2.2. Discussion of possible biases 

It is very important to evaluate the possible sources of biases, 
which may affect the spectral classification, coming from the 
use of various diagnostics at different redshift ranges. 

As the [Oll]A3727/H)3 line ratio is less accurate than 
[Nii]A6584/Ha, or [Sii]AA6717H-6731/Ha, to distinguish be- 
tween star-forming and Seyfert 2 galaxies, the blue diagnostic is 
defined with an error domain (see the dashed curves in Fig. |4]i, 
inside which individual galaxies cannot be safely classified. We 
thus introduce two new categories of galaxies; candidate star- 
forming galaxies and candidate AGNs, which fall respectively 
in the lower-half or the upper-half part of the error domain. 

We know however that AGNs are in minority in the universe. 
Thus we emphasize that /) still no reliable classification can be 
performed for any individual galaxy falling inside the error do- 
main of the blue diagnostic, ;/) the search for AGNs is highly 
contaminated in the candidate AGN region, but Hi) any study 
involving statistically significant samples of star-forming galax- 
ies should not be biased when candidate star-forming galaxies 
and candidate AGNs are included. As it will be shown in Sect.|5] 
and|6] this ends up to a negligible bias on the derived metallici- 
ties. 

As for the blue diagnostic, we have defined an error do- 
main for the Hj3 diagnostic in the following range: 0.4 < 
log([Oiii]A5007/Hj3) < 0.6. Star-forming galaxies or AGNs 
falling inside this domain are classified as candidates. 

One can see in Table. [3]that the red diagnostic is the only one 
able to find LINERs. Indeed, LINERs fall in the star-forming 
galaxies region with the blue or Hj3 diagnostic, while they fall in 
the Seyfert 2 region with the Ha diagnostic. We know however 



Another bias could come from the population of composites 
galaxies, i.e. for which the ionized gas is produced by both an 
AGN and some s tar-forming regions. Whe n looking for AGNs 
in the SDSS data. iKauffmann et alj (l2003al) have defined a new, 
less restrictive, empirical separation between star-forming galax- 
ies and AGNs (see the dashed curve in Fig. O. They have then 
classified as co mposites all ga l axies between this new limit and 
the old one bv iKewlev et al. ( 20011). This res ult is confirmed 
by theoretical modeUng: IStasiiiska et alJ (l2006h have found that 
composite g alaxies are indeed fallin g in the region between 
the curves of'K auffmann et a l.' (2003a') and lKewlev et al.l (1200 ll) 
in the [Oiii]A5007/Hj3 vs. [NiijA6584/Ha diagram. Moreover, 
they have found that composite galaxies fall in the star- 
forming galaxy region in the other diagrams ([Olll]A5007/Hj3 
vs. [Sii]AA6717H-6731/Ha or [On]A3727/H/3). 



It is difficult to evaluate the actual bias due to contamina- 
tion by composite galaxies, in the red or blue diagnostics that 
we use in this study. This difficulty comes mainly from the fact 
that, if comp osite galaxies actu a lly fall in the composite region 
as defined by IKauffmann et al.l (l2003ah . not all galaxies inside 
this region are necessarily composites. A large majority of them 
might be normal star-forming galaxies. One way to evaluate the 
contamination by composite galaxies is to look for star-forming 
galaxies with an X-ray detection. Such work has been performed 
with zCOSMOS data with similar selection criteria than the wide 
sample: Bongiorno et al. (in preparation) have found that the 
contamination of star-forming galaxies by composites is approx- 
imately 10%. 



3. Estimation of the stellar masses and absolute 
magnitudes 

3.1. The Bayesian approach 

The stellar masses are estimated by comparing the observed 
Spectral Energy Distribution (hereafter SED), and two spectral 
features (H5w absorption line and D„ (4000) break), to a library 
of stellar population models. The use of two spectral features 
reduces the well-known age-dust-metallicity degeneracy in de- 
termining the mass-to-light ratio of a galaxy. Compared to pure 
photometry, Balmer absorption lines are indeed less sensitive to 
dust, while the D„(4000) break is less sensitive to metallicity, 
both being sensitive to the age. 

One observation is defined by a set Fi of observed fluxes in 
all photometric bands, a set aFi of associated errors, a set /,/ 
of observed indices, and a set (T/,/of associated errors. The 
of each model, described by a set Fj^ of theoretical photometric 
points, and a set of theoretical indices, is calculated as: 



X 



aF? 



(3) 



where A is the normalization constant that minimizes the x^. 
Note that the normalization constant has not to be applied to the 
spectral indices as they are already absolute. 

Still for one observation, the set of xj calculated on all mod- 
els are summarized in a PDF (Probability Distribution Function) 
which gives the probability of each stellar mass, given the under- 
lying library of models (the prior). Each stellar mass M^, is as- 
signed a probability described by the following normalized sum: 



P{MA{Fi}) 



L,- 5(M, - A ■ Mi) ■ exp (-^7/2) 



(4) 



where Mi is the stellar mass of the model associated to j?. In 
practice, the PDF is discretized into bins of stellar masses. Our 
stellar mass estimate is given by the median of the PDF . 

This method, based on the Bayesi an approach, has been in- 
troduced by the SDSS collaboration ("Kauffm ann et al.l l2003bl: 
[Tremonti et al. 2004; Brinchmann et al. 2004) in order to carry 
SED fitting estimates of the physical properties of galaxies, and 
starts now to be widely used on this subject. Its main advantage 
is that it allows us, for each parameter, to get a reliable estimate, 
independently for each parameter, which takes all possible solu- 
tions into account, not only the best-fit. Thus, this method takes 
into account degeneracies between observed properties in a self- 
consistent way. It also provides an error estimate of the derived 
parameter from the half-width of the PDF . 

3.2. Description of the models 

We use a library of theoretical spectra based on 
(iBruzual & Charlod l2003l hereafter BC03) stellar popula- 
tion synthesis models, calculated for various star formation 
histories. Unless many other studies which only use a standard 
declining exponential star formation history, we have used 
an imp roved grid i ncluding also secondary bursts (stochastic 
hbrary, ISaUm et al.l |200 5; Gallazzi et al. 2005). The result of 
the secondary bursts, compared to previous methods, is mainly 
to get higher masses as we are now able to better reproduce 
the colors of galaxies containing both old and recent stellar 
populations. When one uses a prior with a smooth star formation 
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Figure 6. Comparison of the stellar masses of star-forming 
galaxies in our sample, obtained with BC03 and CB07 models, 
using both photometry and spectroscopy. The latter provides a 
better treatment of the TP-AGB stars in the near-infrared. The 
top panel shows the difference between the two masses, the bot- 
tom panels show the associated histogram (left) and the his- 
togram of the differences between the two error estimates (right). 
The slope of the difference (s) is shown in the plot. The mean 
value (m) and the rms (r) of the two histograms are given in the 
plots. 



history, such objects are better fitted by young models which 
correctly reproduce the recent star formation history and the 
colors in the bluer bands. But such models fail to reproduce the 
underlying old stellar population which affects the colors of the 
redder bands, therefore they lead to an underestimate of the final 
stellar masses by a factor « 1 .4 (see Pozzetti et al. 2007, for a 
discussion of the different methods). 

Masse s derived with t his method are also relative to the cho- 
sen IMF (IChabriedl2063l) . The models include self-consistent 
two-component dust corrections (ICharlot & Fallll2000h . 

In addition to the secondary bursts, the models used in 
this s tudy also include a n ew treatment of the TP-AGB 
stars dMarigo & Girardil |2007|) (Chariot &. Bruzual, in prepa- 
ration, hereafter CB07), in order to provide a better fit of 
the near-infrared pho tometric bands (see also iMarastonI 120051 
iMaraston et al.ll2006l) . We have derived the stellar masses of our 
whole sample of galaxies with BC03 and CB07 models. Figure|6] 
shows the results of this comparison. As expected, the masses 
derived with CB07 models are smaller, since the near-infrared 
flux of TP-AGB stars was underestimated in BC03 models. The 
mean shift is —0.07 dex, with a dispersion of 0.11 dex, which 
corresponds to an overestimate of the stellar masses by a factor 
1.2 with BC03 models. Their is also a small trend for more 
massive galaxies to be more affected by the difference between 
the two models. This compa rison has been studied in more de- 
tails by lEminian et al.1 (l2008l) . 

The same comparison have been performed on masses de- 
rived only using photometry: the mean shift is —0.09 dex, with a 



dispersion of 0.19 dex. The spectral indices are not expected to 
vary much with the new treatment of TP-AGB stars introduced 
by CB07 models. It is thus expected that masses computed only 
with photometry vary more from BC03 to CB07 models than 
those computed also with spectral indices. 

We note that the stellar masses derived with BC03 models 
and a smooth star formation history are very similar to the ones 
derived with CB07 models and secondary bursts. Stronger ef- 
fects have to be expected on derived star formation rates or mean 
stellar ages. 

Figure |7] now compares the stellar masses obtained with 
CB07 models when using only photometry, only spectroscopy, 
and both photometry and spectroscopy. Figure |7Ja) shows the 
comparison between the mass obtained using only photometry, 
and the one obtained by adding the two spectral indices men- 
tioned above. The derived mass and associated error are not sig- 
nificantly biased towards the two different methods (photometric 
masses are —0.05 dex lower). Nevertheless, we have compared 
the values in both cases and found that it never increases 
when we add the spectral indices in the analysis. Moreover, the 

obtained in this latter case is lowered by a factor whose mean 
value is 2.3. The dispersion of 0.15 dex in the comparison be- 
tween the two masses has thus to be taken into account as the 
minimum uncertainty of our stellar masses. 

Figure |7Jb) shows what would happen if only spectral in- 
dices were used. We see that the derived mass does not suffer 
from any bias, but is subject to a large dispersion of 0.39 dex (a 
factor 2.45). A slope of —0.25 dex/decade is found, but is not 
significant given the dispersion. The comparison of the derived 
errors shows that this additional uncertainty is mostly taken into 
account: both measurements have to account for a minimum un- 
certainty of 0.15 dex and the measurement obtained using only 
spectroscopy show a 0.13 dex additional uncertainty. 

The reason of this comparison is that lTremonti et all (l2004l) 
have not used full SEDs to derive stellar mass-to-light but only 
two spectral indices: D„(4000) and H5a, the mass being then 
scaled from z-band luminosity. We note that their spectral in- 
dices have been measured at a much better signal-to-noise and 
resolution on SDSS spectra as compared to V VP S spectra. We 
conclude from Fig.lTJb) that lTremonti et all (l2004l) masses com- 
puted only with spectroscopy are directly comparable to our final 
masses computed also with photometry, provided that we take 
scale them by —0.07 dex and —0.056 dex to account respec- 
tively fo r the difference bet ween BC03 and CB07 models, and 
betweenl&ouii (|2001|) and lChabriej (|2003|) IMF. 

3.3. Defining volume-limited samples 

Figure [8] shows the stellar masses obtained, as a function of red- 
shift, for the galaxies of our various samples (wide or deep, 
all galaxies or star-forming galaxies only). As expected in a 
magnitude-selected sample, the minimum stellar mass which is 
detected depends on the redshift, because of the decrease of 
apparent magnitude when the distance increases. This effect is 
well-known as the Malmquist bias. The upper-limit of the mass 
distribution is in part due to the bright cut-off in the magnitude 
selection, while in addition to that a larger volume is observed at 
higher redshift which allows a better sampling of the upper mass 
region. The VVDS is a constant solid-angle survey: the observed 
volume increases with distance and rare objects, like very mas- 
sive galaxies, get more chances to be detected. 

There are two ways to account for these two biases. The first 
one is to use the Vmax technique: each object is weighted by the 
inverse of the maximum volume in which it could be observed. 



given the selection function of the survey. The result is a den- 
sity which is statistically corrected for the selection effects, but 
may suffer from evolution effects. Another, simpler, technique is 
to build volume-limited samples: a redshift range is defined and 
only objects which may be observed in the entire volume be- 
tween these two redshifts (given the selection function) are kept. 
We use this last technique in our study. 

In practice, for magnitude-selected surveys, one has to cal- 
culate the minimum stellar mass which is detected at the upper 
limit of the redshift range, and then to throw out the objects be- 
low this limit since they are not observed in the entire volume. 
Nevertheless the minimum mass which may be detected at a 
given redshift, for a given limiting apparent magnitude, depends 
on the stellar mass-to-light ratio. One way to solve this issue is 
to calculate, for each galaxy, the mass Mj*^^^ it would have if its 
apparent magnitude / was equal to the selection magnitude /sei 
(i.e. 24 for the deep sample, and 22.5 for the wide sample): 

A^llm=A^*+0.4-(/-/sel) (5) 

The result is a distribution of minimum stellar masses which 
reflects the distribution of stellar mass-to-light ratios in the sam- 
ple. The red, green, and blue curves in Fig. [8] show respectively 
the 50%, 80% and 95% levels for the cumulative sum of this 
distribution. The limiting masses derived for the 4 different sam- 
ples mentioned above are given in Table 2] As expected, we see 
that the limiting mass increases with redshift and is lower for 
the deep sample. The derived completeness leve ls on the deep 
sampl e are similar but higher to the ones found by lMeneux et al.l 
(I2008h on the same data. The difference is due to the different 
cosmology used in this paper, i.e. h = I. 

We see also that the limiting mass is lower when we consider 
only star-forming galaxies, compared to the whole samples. This 
is due to star-forming galaxies having lower mass-to-light ratios. 
This shows that if a sub-sample has a smaller range in mass-to- 
light ratios, then a less conservative mass limit can be adopted 
than if a single fixed maximum mass-to-light ratio is taken. 

Finally, we note that the wide sample contains more massive 
objects since it covers a larger volume of universe. 

3.4. Absolute magnitudes 

As a by-product of the estimation of stellar masses, we have also 
calculated the fc-corrected absolute magnitudes, in the rest-frame 
Z?-band, using the Bayesian approach and the CB07 models. The 
results for star-forming galaxies are shown in Fig.|9]as a function 
of redshift. Because of the ^-correction, the maximum absolute 
magnitude which may be detected, as a function of redshift, is 
also affected by the stellar mass-to-light ratio. The absolute mag- 
nitude in rest-frame Z?-band is affected by the ratio between the 
flux in this band, and the flux in the observed /-band, which ra- 
tio depends on the type of the galaxy and on redshift We also 
note that as observed /-band is fairly close to rest-frame B-band, 
the variation in A:-correction is small compared to the variation in 
mass-to-light ratio as can be seen by the reduced spread between 
the different completeness level curves in Fig.[8]as compared to 
Fig.|9] 

The limiting B-band absolute magnitudes, associated to var- 
ious redshift ranges and ^-correction completeness levels, are 
given in Table |5] We note that the absolute magnitudes used 
in this study may d iffer from previou s works performed on the 
same dataset (e.g. Illbert et ani2005h . For self-consistency, we 
have chosen to calculate the absolute magnitudes with the same 
models than the ones used for the stellar masses, namely CB07 
models with secondary bursts. 



0.5 



s = 0.03 




250 
200 
150 
100 
50 




1 1 1 1 1 1 1 1 1 


1 L 


- m = 0.05 




; r = 0.15 


1 : 











10 

Mass (Phot) 



I- 200 

e 

^ 100 



11 



12 



-1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1. 


r m=-0.00 




- r=0.06 














-1 -0.5 0.5 1 
A(Mass) (tinal-Phot) 



-0.4 -0.2 0.2 0.4 
A(aMass) (final-Phot) 



(a) 




40 
30 

0) 

E 20 



m=-0.00 
- r = 0.39 





10 11 
Mass (final) 



60 



20 



J 1 1 1 1 1 1 1 1 

- m = 0.13 






1 1 1 1 i_ 


- r = 0.10 






























-1 -0.5 0.5 1 
A(Mass) (Spec-final) 



-0.4 -0.2 0.2 0.4 
A((7Mass) (Spec -final) 



(b) 



Figure 7. Same legend as in Fig.|6] except that we now compare the stellar masses that we obtain with CB07 models with photometry 
and spectroscopy (final) to the ones obtain; (a) using only photometry, (b) using only spectroscopy. Only objects with a sufficient 
number of photometric bands and spectral indices (i.e. more than two of each) are used in these plots. 



Table 4. Limiting masses (in logarithm of solar masses) of our various samples in order to define volume-limited samples, as a 
function of the redshift ranges and for various mass-to-light completeness levels. 
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Table 5. Limiting B-band absolute magnitudes (AB system) of our wide and deep samples of star-forming galaxies, in order to 
define volume-limited samples, as a function of the redshift ranges and for various A:-correction completeness levels. 



z<0.5 z<0.7 z<0.9 

50% 95% 50% 95% 50% 95% 
wide -19.1 -19.6 -19.9 -20.3 -20.5 -20.9 
deep -17.5 -18.3 -18.2 -18.9 -18.6 -19.3 



4. Estimation of the metallicities 

In this study, the metallicities are estimated as the gas-phase oxy- 
gen abundances 12 + log (O/H). 

4.1. The empirical approach 

Strong-line empirical cahbrators may be used to compute the 
gas-phase metallicities from the available measured emission 
lines, e.g. Ha, Hjg, [OiiiU5007, [N iiU6584, [Oii]A3727. As 
shown by iLamareille et al.l (l2006bh . the [Nll]A6584 and Ha 
emission lines are reliably deblended for sufficiently high signal- 
to-noise ratio even at the resolution of the VVDS. 

We use the N2 c alibrator up to z 0.2 (iDenicolo et alj2002t 
Ivan Zee et al ] ll998h . the 03N2 calibrator for the redshift range 



0.2 < z < 0.5 (iPettini & Pagelll2004. and the Ri t, calibrator for 
the redshift range 0.5 < z < 0.9 (lMcGaughlll99ll) . 

At high redshifts, where nitrogen or sulfur lines are not ob- 
served, the metallicity is degenerate because of radiative cool- 
ing at high metalhcities which lowers the intensities of oxygen 
lines. Thus, low oxygen-to-hydrogen line ratio can be interpreted 
as either low oxygen abundance or high oxygen cooling, i.e. 
high oxygen abundance. The degeneracy of the R23 calibrator is 
broken using the B-band absolute magnitude and the refere nce 
luminosity-metallicity relation (see ILamareille et^l2006aL for 
a detailed discussion of this method). 

As clearly shown bv lKewlev & ElUsonl (|2008|) . the difficulty 
of using different metallicity calibrators at different redshifts 
comes from the huge differences, up to ±0.7 dex, which exist 
between them. Therefore, we have decided to renorm alize the 
different calibrators to the ICharlot & Longhettil (1200 iL hereafter 
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Figures. Stellar mass estimates (in logarithm of solar masses), as a function of redshift, of the galaxies in our samples; left; the 
wide sample, right; the deep sample, top; all galaxies (except broad-line AGNs), bottom; only star-forming galaxies (see Sect. 12.2b . 
The red, green and blue curves illustrate respectively the 50%, 80% and 95% mass-to-light completeness levels (see text for details). 



CLOl) calibrator, which has been u sed to estimate the mass- 
metallicity relation from SDSS data (iTremonti et al.ll2004. The 
correc tion formulas are found in Table B3 of iKewlev & EllisonI 
(120081) . 

Finally, a consequence of the small spectral coverage is also 
that Ha and Hj3 lines are not observed together for most red- 
shifts, making difficult the estimation of the dust attenuation 
fr om the standard B aimer decrem ent r nethod. As already show n 
by iKobulnicky & PhilUpsI (l2003l) . and lLamareille et al.1 il2006al) . 
metallicities can be reliably derived using equivalent width mea- 
surements instead of dust-corrected fluxes. 



Nevertheless, iLiang et aP (l2007h have shown that the 
[Olll]A5007/[Oll]A3727 emission line ratio, which is used in 
R2i calibrator, does not give the same value if it is calculated us- 
ing equivalent widths or dust-corrected fluxes; the ratio between 
the two results depends on the differential dust attenuation be- 
tween stars and gas, and on the slope of the stellar continuum. 
They have found that the derived metallicities are affected by a 
factor ranging from —0.2 to 0.1 dex, but they also showed that 
this factor has a mean value of only —0.041 dex. It is thus not 
significant for the estimation of the mass-metallicity relation on 
statistical samples. Nevertheless, in order to avoid any possible 
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Figure 9. Absolute magnitudes estimates in the rest-frame B-band (AB system), as a function of redshift, of the star-forming galaxies 
of the wide (left) and deep (right) samples. The red, green and blue lines illustrate respectively the 50%, 80% and 95% ^-correction 
completeness levels (see text for details). 



bias on a different sample, we decided to use their corrective 
formula based on the D„ (4000) index. 

4.2. The Bayesian approach 

The metallicities may also be estimated using again the Bayesian 
approach. The relative fluxes of all measured emission lines 
are compared to a set of photoionization models, which pre- 
dicts the theoretical flux ratios given four parameters: the gas- 
phase metallicity, the ionization level, the dust- to-metal ra- 
tio and the reddening dCharlot & Longhettill200Th . The CLOl 
mod els are based on popul ation synthesis for the ionizing 
flux dBruzual & Charloj(2003h . emission line modeling (Cloudy, 
iFerland 20011) and a two-component dust attenuation law 
tCharlot & Fallll2000l) . 

We calculate the of each model and summarize them in 
the PDF of the metallicity using a similar method than described 
above in Eq.[3]and|4] but now applied to emission-line fluxes in- 
stead of photometric points. Only the emission lines with enough 
signal-to-noise (S/N>4) are used in the fit. This method is ap- 
plied directly on observed line fluxes: the correction for dust at- 
tenuation is included self-consistently in the models. 

The O/H degeneracy at high redshift ends up with double- 
peaked PDFs: one peak for the low metallicity solution and an- 
other one at high metallicity. However, other information such 
as dust extinction, ionization level, or star formation rate help 
affecting to the two peaks two different probabilities. Thus, we 
solve the degeneracy by fitting two peaks in the PDFs, and by 
keeping t he one with the highest probability. As already dis- 
cussed in iLamar eille et al.l (l2006ah . this method cannot be used 
to chose the metallicity of one single galaxy, but it can be used 
statistically to derive a mean metallicity, e.g. as a function of 
mass. 

In this study, given the relatively low spectral resolution and 
spectral coverage of our spectra, we use the CLOl method only 



as a check of the quality of empirically-calibrated metallicities. 
The three main advantages of checking our results with the CLOl 
method are: ( i) the use of one unique calibrator for all redshift 
ranges; (ii) the different method to deal with dust extinction, 
i.e. the self-consistent correction instead of the use of equivalent 
widths; (Hi) the different method to break the O/H degeneracy, 
i.e. the fit of the double-peaked PDF instead of the use of the 
luminosity. 

Fig. [To] shows the comparison between the empirically- 
calibrated and the CLOl metallicities. Apart from a small num- 
ber of outliers (7%), the two results are in good agreement within 
a dispersion of approximately 0.22 dex. This dispersion can be 
compared with the intrinsic dispersion of the mass-metallicity 
relation, which is also of ~ 0.22 dex. We have also calculated 
the mean of the error bars for the two methods: 0.06 dex for the 
empirical method, and 0.18 dex for the CLOl method. 

Let us summarize all the contributions to the global dis- 
persion, (i) The empirically-calibrated errors have been esti- 
mated from basic propagation of the line measurement errors. 
We know from tests performed on duplicated observations that 
platefit_vimos errors are not underestimated. They reflect the 
negligible contribution of noise to the global dispersion, (ii) 
Conversely, the CLOl errors have been estimated from the width 
of the PDF. They also reflect the additional dispersion due to the 
degeneracies in the models, which are taken into account thanks 
to the Bayesian approach. (Hi) Finally the difference between 
the mean CLOl error and the global dispersion reflects real vari- 
ations in the physical parameters of galaxies, which are not taken 
into account in both models. 

The comparison between the two methods have shown 
that: ( i) the various empirical calibrators used at different red- 
shift ranges give consistent results with the CLOl metallicities, 
thanks to the correction formulas provided by lKewley & EllisonI 
( |2008|) ; (ii) the empirical correction for dust attenuation us- 
ing equivalent widths, and the correction formula provided by 
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Figure 10. Comparison between empirically-calibrated (x-axis), and CLOl metallicities (y-axis) for our whole sample of star- 
forming galaxies. The red solid line shows the y ~ x relation, while the blue dashed line is a least-square fit to the data points. 
The parameters of the fit are given in the two plots (a: slope, b: zero-point). 



iLiang et alJ (l2007h . is consistent with CLOl's results, obtained 
directly from observed line fluxes and a self-consistent redden- 
ing correction; ( Hi) apart from a relatively small number of out- 
liers, the O/H degeneracy is correctly broken using the luminos- 
ity diagnostic with the empirical method. 



4.3. Possible biases 

Fig.[TO]also indirectly shows that the metallicities obtained with 
N2 and 7?23calibrators at different redshifts are directly compara- 
ble, even if there is no galaxy in our sample where these two cal- 
ibrators can be applied at the same time. Indeed we see that both 
N2 and /?23metallicities are in agreement with CLOl metallici- 
ties, apart from the small number of outliers. Moreover we have 
investigated the possible biases introduced by the use of differ- 
ent sets of lines, at different redshifts, with the CLOl method: we 
have confirmed that this method gives, as expected, comparable 
results at different redshifts. 

Another possible bias comes from the fact that the slit used 
for the spectroscopy not necessarily covers all the light of the ob- 
served galaxy. We have taken the spatial extension of the galax- 
ies calculated from photometry, and compared it to the width 
of the slit, i.e. one arcsecond. This gives the aperture fraction 
of the galaxy, expressed in percent. Fig. fTTl bottom) shows the 
distribution of the aperture fractions for the star-forming galax- 
ies in our sample (not that this information was not available 
for the CDFS sample). The mean aperture fraction is 52%, and 
we see that the large majority of objects has an aperture frac- 
tion grea ter than 20%, which is the minimum aperture fraction 
given by iKewley. Jansen. & Gelled (l2005h in order to minimize 
aperture effects on derived parameters hke the metallicity or the 
star-formation rate. 
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Figure 11. Bottom: distribution of the aperture fractions (in per- 
cent) for the star-forming galaxies in our sample. The dashed 
histogram shows the distribution for the galaxies at low redshift 
(z < 0.5) with a high metallicity (12 + log(0/H) > 8.8). Top: 
evolution of the apertures fraction as a function of redshift. The 
solid line is the fit to the data points. 



In Fig.[TTl we also observe the expected correlation between 
redshift and the aperture fraction, both increasing at the same 
time. The aperture fractions show an increase of 35% per unit- 
redshift. Since the aperture fractions are already high enough 



at low redshift not to affect the derived metallicities, and since 
they increase with redshift, this slope is likely to have a marginal 
effect on the derived evolution of metallicity as a function of 
redshift. 



5. The luminosity-metallicity relation 

5.1. Study of the derived fits 

We study in this section the relation between the rest-frame B- 
band luminosity and gas-phase metallicity, for the star-forming 
galaxies of our wide and deep samples. Here we use the metal- 
Ucities estimated with the empirical approach, and renormal- 
ized to the CLOl method (see Sect. 14. it . Figure [T2l shows the 
luminosity-metallicity relation of the wide and deep samples in 
three redshift ranges: 0.0 < z < 0.5, 0.5 < z < 0.7 and 0.7 < z < 
0.9. The results are compared to th e luminosity-metallicity rela- 
tion in the local universe derived by lLamareille et all (l2004l) with 
2dFGRS data, and renormalized to the CLOl method. In order to 
do this comparison, we have done a linear fit to the data points. 

As already shown before, the luminosity-metallicity relation 
is characterized by an non-negligible dispersion of the order of 
« 0.25 dex (higher than the one of the mass-metallicity relation). 
Therefore, the method used to perform the fit has a huge impact 
on the results, and it is mandatory to use the same method before 
doing comparisons between two studies. Th us, we have used the 
same method than Lamareille et al. I (12004 . i.e. the ols bisector 
fit dlsobe et al ] ll990ir starting from the 50% ^-correction com- 
pleteness level (see Sect. 13.4b . The results are shown in Fig. [12] 
and in Table |6] 

The slope is by ~ 2(7 steeper than the one of the reference re- 
lation in the local universe, this slope being similar in all redshift 
ranges (« —0.75), except 0.5 < z < 0-7 in the deep sample. As 
already shown by iLamareille et al.l (|2004|) who have compared 
the slopes of the luminosity-metallicity relation for low and high 
metallicity objects, this slope is steeper f or high metallicit y ob- 
jects. This effect has been confirmed by iLee e t al.^ ( 200 6|) who 
have e xtended the mass-metallicity relation of Tremonti et al.l 
(12004 to lower masses and metallicities and found a flatter 
slope. We emphasize that we are not discussin g here the sat- 
uration at 12 + log(0/H) « 9.2 as observed by iTremonti et alJ 
(12004 . and which might be confused with a flatter slope at very 
high metallicity. 

Thus, given that we expect a steeper slope at high metallicity, 
we analyze our results as a lack of data points in the low metal- 
licity region: the wide sample does not go deep enough to detect 
low-metallicity and low-luminosity objects above the complete- 
ness limit. Moreover, in the 0.0 < z < 0.5 redshift range of wide 
and deep samples, the S/N>4 cut on emission lines introduces 
a bias towards high metallicity objects: the faint and blended 
[Nll]A6584 line, used to compute metallicities in this redshift 
range, becomes rapidly undetectable at low metallicities. 

The apparently high number of galaxies with low redshift, 
low luminosity and high metallicity seen in Fig. fT2l left) may 
be explained by two effects. First, as said in previous paragraph, 
the [Nll]A6584 line is more likely to be observed for high metal- 
licity objects. There is consequently a lack of objects with low 
luminosity and low metallicity which makes the distribution ap- 
parently wrong (as compared to previous studies). Second, these 
objects at low redshift may have lower aperture fractions, which 
may introduce a bias towards higher metallicities given that we 
preferentially observe their central part. Indeed, we have plot 
the distribution of the aperture fractions for objects with a low 



redshift and a high metallicity in Fig. fTTT bottom. dashed his- 
togram): they are clearly lower, with a mean of 40%. 

We recognize that the results obtained by doing a direct lin- 
ear fit to data points suffer from a drawback: the correlations 
between metallicity and luminosity shown in Fig. [T2] are weak. 
This weakness is characterized by the non-negligible error bars 
for the slopes and the zero-points of the relations provided in 
Table|6l and also by low Spearman correlation rank coefficients: 
of the order of —0.3 for the bottom-center panel, and— 0.1 for 
the other panels. 



5.2. Global metallicity evolution 

Despite the weakness of the correlations found in our data, we 
know from previous studies that the luminosity-metallicity ex- 
ists. Thus we can use the existence of this relation as an assump- 
tion and find new results. 

In order to quantify the evolution of the luminosity-relation, 
we can also derive the mean evolution of metallicity. This is done 
with the additional assumption that the slope of the luminosity- 
relation remains constant slope at zero-order The results are 
shown in Fig. [12] and in Table [6] As expected, the evolution 
is barely significant in the 0.0 < z < 0.5 redshift range. In the 
0.5 < z < 0.7 and 0.7 < z < 0.9 redshift ranges, the evolution is 
str onger. It is simila r in the wide sample to what has been found 
by lLamareille et al.l (l2006al) . In the 0.5 < z < 0.7 redshift range 
of the deep sample, the results are similar with both methods, 
which confirms the hypothesis that a steeper slope is found only 
when low-metallicity points are not included in the fit. 

Finally, the comparison of the wide and deep samples show 
a stronger evolution of the metallicity with redshift in the wide 
sample. This result tells us that the slope does not actually 
remain constant, and that the evolution of the metallicity is 
stronger in more luminous objects. At z 0.76, galaxies with an 
absolute B-band magnitude of ~ —20.1 have —0.32 dex lower 
metallicities than galaxies of similar luminosities in the local 
universe, while galaxies with an absolute B-band magnitude of 
^ —21.2 have —0.65 dex lower metallicities. 

We have checked that the results are stable when break- 
ing the O/H degeneracy at high redshift with a different refer- 
ence relation (e.g. lower ed in metalUcity). As stated before by 



ILamareille et al.l72006al) . the use of a lower reference relation to 
break the degeneracy only changes the metallicities of the galax- 
ies in the intermediate region (12 + log(0/H) « 8.3), thus non- 
significantly changing the whole luminosity-metallicity relation. 

We have also checked the effect of not including the candi- 
date star-forming galaxies in the fit. As stated in Sect. 12.2! the 
contamination of candidate star-forming galaxies by AGNs is 
less than 1%, and is thus not expected to affect significantly the 
luminosity-metallicity relation. Therefore, the results of not in- 
cluding the candidate star-forming galaxies are only given as in- 
formation in Table [6] We have also checked the effect of includ- 
ing candidate AGNs. In both cases, the effect on the luminosity- 
metallicity relation is very small. 

We need to disentangle between the two effects which make 
these objects more luminous, i.e. a higher mass or a lower mass- 
to-light ratio, which one is responsible of the stronger evolution 
of the metallicity. Thus, we now study the mass-metallicity rela- 
tion in the next section. 




Figure 12. Rest-frame B-band luminosity-metallicity relation for the wide (top) and deep (bottom) samples, for three redshift ranges: 
from left to right 0.0 < z < 0.5, 0.5 < z< 0.7 and 0.7 < z < 0.9. The metallicities have been estimated using the empirical app roach 
(see Sect. 14. it . The solid line shows the luminosity-metallicity relation at low redshift derived bv iLamareille et al] (|2004|) . and 
renormalized to the CLOl method. The short-dashed line shows the fit to the data points using the ols bisector method (see text). 
The long-dashed lines shows the fit to the data points assuming a constant slope. 

Ta ble 6. Evolution of the luminosity-metallicity relation for the wide and deep samples. The reference relation is the one obtained 
by ILamareille et aP (l2004l) with 2dFGRS data and renormalized to the CLOl method. In each redshift ranges, we give the results of 
the ols bisector fit (slope and zero-point) and of the constant-slope fit (mean shift, see Fig.[T2]i. together with the mean redshift and 
magnitude, and the dispersion of the relation. The metallicity shift is given in three cases: a) using only star-forming galaxies, b) 
adding candidate star-forming galaxies, and c) adding also candidate AGNs. 
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6. The mass-metallicity relation 

6.1. Global metallicity evolution 

We now derive the mass-metallicity relation for the star-forming 
galaxies of the wide and deep samples, using the empirical ap- 
proach for computing the metallicities, and in three redshift 
ranges: 0.0 < z < 0.5, 0.5 < z < 0.7 and 0.7 < z < 0.9. The re- 
sults are shown in Fig.[T3]and Table|2l and are compared to the 
re ference mass-metallici ty relation in the local universe, derived 
by lTremonti et alJ (l2004l) with SDSS data. The reference relation 



has been shifted in mass, in order to take into account the effect 
of using different models (see Sect. 13.21 ). 

As in Fig.[l2] the data points shown in Fig. [13] do not show 
strong Spearman rank correlation coefficients. We thus skip the 
step of doing a fit to these data points. We directly probe the 
global evolution in metallicity of star-forming galaxies com- 
pared to the reference relation, doing the likely assumption that 
the mass-metallicity relation exists also at high redshift. To do 
so, we calculate the mean shift in meta llicity by fitting to th e 
data points the same curve than Eq. 3 of iTremonti et al.l (l2004l) . 
allowing only a different zero-point. The zero-order assumption 
is indeed that the shape of the mass-metallicity relation does 




Figure 13. The mass-metallicity relation of star-forming galaxies for the wide (top) and deep (bottom) samples, for three redshift 
ranges: from left to right 0.0 < z < 0.5, 0.5 < z < 0.7 and 0.7 < z < 0.9. The metallicities have been e stimated using the emp irical 
approach (see Sect. 14. lb . The solid curve shows the mass-metallicity relation at low redshift derived bv lTremonti et al.l (|2004|) . The 
long-dashed curves shows the fit to the data points assuming that the SDSS curve is only shifted down in metallicity. 



Table 7. Evolution of the mass-metallicity relation for the wide 
and deep samp l es. Th e reference relation is the one obtained by 
iTremonti etal] (|2004 with SDSS data and the CLOl method, 
renormalized in mass. In each redshift ranges, we give the mean 
redshift, stellar mass, and metallicity shift (assuming that the 
shape of the relation remains constant, see Fig. [13]), and the dis- 
persion of the relation. The metallicity shift is given in three 
cases: a) using only star-forming galaxies, b) adding candidate 
star-forming galaxies, and c) adding also candidate AGNs. 
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not vary with redshift. The results are shown in Fig. [T3] and in 
Table [T] The fit is performed above the 50% mass-to-light com- 
pleteness level. 

There are some caveats in the interpretation of the mass 
and metallicity evolution of the deep sample, because of selec- 
tion and statistical effects. The mean observed stellar seams to 
decrease with redshift, which is in contradiction of what one 
would expect from the Malmquist bias. First, we find higher 
stellar masses than expected in the lowest redshift bin. The sam- 



ple is actually not complete down to the limiting mass: lower 
mass galaxies would have lower metallicities, and low metallic- 
ities are indeed difficult to measure because of the blending of 
[Nll]A6584 and Ha lines. Conversely, the highest redshift bin 
has a lower mean mass than expected, which is more probably 
due to a statistical effect because of the small solid angle of the 
deep sample. This later effect also explains why the metallicity 
evolution seems smaller in the highest redshift bin. 

As observed on the luminosity-metallicity relation in previ- 
ous section, we clearly see a stronger metallicity evolution in 
the wide sample than in the deep sample. The wide and deep 
samples span interestingly different ranges in masses, but are 
otherwise identical. This effect thus shows that the most mas- 
sive galaxies have experienced the most significant evolution 
in metallicity. At z 0.77, galaxies at 10^ '* solar masses have 
—0.18 dex lower metallicities than galaxies of similar masses 
in the local universe, while galaxies at lO^"'^ solar masses have 
—0.28 dex lower metallicities. We therefore conclude that the 
shape of the mass-metallicity relation varies with redshift, so that 
it was flatter in an earlier universe. 

We also remark that like for the luminosity-metallicity re- 
lation, the potential influence of candidate AGNs or candidate 
star-forming galaxies is negligible (see Table|7]i. 

6.2. Evolution of the shape of the mass-metallicity relation 

We now evaluate the evolution of the shape of the mass- 
metallicity relation as a function of redshift, by co-adding a num- 
ber of data points together in bins of mass, thus increasing the 
signal-to-noise ratio. Figure [14] shows the results in three red- 



Figure 14. The mass-metallicity relation of star-forming galaxies for three redshift ranges: from left to right 0.2 < z < 0.5, 0.5 < 
z < 0.6 and 0.6 < z < 0.8. The metallicities have been es timated using the empi rical approach (see Sect. 14. lb . The solid curve shows 
the mass-metallicity relation at low redshift derived by iTremonti et al.l (|2004|) . We show the mean metallicities by bins of stellar 
masses. The blue error bars represent the uncertainty on the mean, while the green error bars represent the dispersion of the data 
points. The 50% mass-to-light completeness levels for the deep and wide samples (see Table|4ll are shown as vertical dotted Unes. 




Figure 15. The mass-metallicity relation of star-forming galax- 
ies for two redshift ranges: 0.5 < z < 0.6 (dark grey) and 
0.6 < z < 0.8 (light grey). The metallicities have been estimated 
using the empirical approach (see Sect. 14.1b . The solid curve 
shows the mass-meta llicity relation at low redshift derived by 
ITremonti et alJ (l2004l) . This relation has also been linearized in 
the range 8.9 < lo^j Mi^/Mp,) < 10.3 . The s hort-dashed curve is 
the relation derived by lSavaglio et al ]('2005) at z ^ 0.7, the long- 
dashed line is the relation derived by fErb et al. (2006) at z - 2.0, 
and the long-dashed- short-dashed line is the relation derived by 
iMaiolino et aP (l2008h at z 3.5. 



shift ranges: 0.2 < z < 0.5, 0.5 < ; 
mean redshifts are respectively z ' 



<0.6and0.6<z<0.8.The 
^ 0.33, z -0.54 and z- 0.67. 



limit to the 50% completeness of the wide sample, above which 
both samples are used. 

The mean masses and metallicities have been estimated in 
respectively 3, 4, and 4 stellar mass bins. This approach is simi- 
lar to the one used at low redshift by |TremontietalJ(l2004 . The 
bins are not equally spaced in mass,but always contain a similar 
number of data points, i.e. — 45 in our case. We have set two 
consecutive mass bins to have 25% of their data points in com- 
mon. We evaluate the errors on the means and the dispersions of 
the data points, which are both shown in Fig. [14] The errors on 
the means are very low thanks to the number of available data 
points. The dispersion in metallicity stays almost constant and 
equal to w 0.22 dex. All results are calculated as the median of 
1000 bootstrap estimates. 

Fig. [TSl shows the evolution of the slope in the two last red- 
shift bins. We confirm that the mass-metallicity relation tends 
to become flatter at higher redshifts: the slope is 0.166 ±0.02 
dex/decade atz ~ 0.54 and 0.150 ±0.01 dex/decade atz ~ 0.67, 
which has to be comp ared to the slope 0.29 dex/decade of the 
ITremonti et alJ (l2004l) curve linearized in the same mass range. 
The slope evolves by —0.25 dex/decade by unit of redshift in the 
range < z < 0.7. These results also confirm the assumption, 
used in previous section, that the mass-metallicity relation exists 
at high redshift. 



6.3. Comparison with previous works 

Fig. [15] shows also the comparison betw een our results, and 
other studies performed at high redshifts dSavaglio et al ] l2005t 
Erb et al. 2006; Maiolino et al. 20081). AU curves have been 
renormalized in stellar masses using the shifts summarize d 
in Table. [8] M etall icities derived bv ISavagUo et alJ (l2005l) . 
lErb et al.] (|2006|) and Maiolino et al. ( 20081) hav e been converted 
respec tively from'Kobulnickv & Kewlev (2004'). rPettini & Pagell 
(|2004 N2 , and .Kewlev &Dopita (2002.) methods to the CLOl 
method. 



The mean masses are respectively 10^'^^, lO^'^"*^ and 10^ ™ solar 
masses. We have joined for this study the wide and deep sam- 
ples: only the deep sample is used from its 50% completeness 



The mass-metallicity relation at redshift z ^ 0.7 is given 
by Eq. 8 of ^avaglio et al. (2005), shifted by 0.47 dex in stel- 
lar mass. The mass-metallic ity relation at redshif t at redshift 
z ^ 2.0 is given by Eq. 3 of ITremonti et alJ (l2004l), shift ed by 
—0.56 dex in metallicity, as found by lErb et alj ( |2006|) . The 
mass-metallicity relation at redshift z — 3.5 is given by Eq. 2 



Table 8. Different setups for different mass-metallicity relations found in the literature, and associated shifts to be applied to their 
stellar masses in order for them to be comparable with our results. This table shows the type of data, the models, the presence or not 
of secondary bursts, and the IMF used to compute the stellar masses. The last column gives the global shift that has to be applied to 
the logarithm of the stellar mass. 
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+0.14 
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*S sta nds for spectroscopy, P stands for p hotometry. 

|T04:lTremonti et alJll2004b: S0 5: Savagl io et al] ( l2005h: E06:lErb et all ( l2006h : M08: lMaiolino et alj ( l2008h : C03: IChabrie3 ( l2003h : K01: lKroupal 
( 120011) : BG03: lBaldry & Glazebrook. t200l) ; S55: ISalpetej a955h . 



of iMaiolino et al.l ( l2008l) . with the parameters given in their 

Table 5. 

Contrary to ISavaglio et al.l (I2005h . we find a flatter slope of 
the mass-metallicity relation at z ~ 1 . Nevertheless, the compar- 
ison of data points shows that their results and ours are in good 
agreement for the highest mass bins. The larger difference comes 
from the lowest mass bins, in which they are probably not com- 
plete. We know indeed that lower mass-to-light ratio galaxies are 
preferentially observed in the lowest incomplete mass bins, and 
that s uch galaxies show smaller mean metallicities (lEUison et alJ 
I2008h . 

The comparison with the data of lErb et all (l2006h . which 
are taken at z 2, is less straightforward: there is a fairly good 
agreement in metallicity with our high-mass end data, but actu- 
ally at a rather different redshift. This could mean that there have 
been very little metallicity evolution from z ~ 2 to z ^ 1, but 
this would be hard to understand when looking also at the other 
results. The m etallicity evolves indeed strongly from z ^ 3.5 
jMaiolino et a l. 2008) , and between z ~ 1 and the local universe 
(our data and ISavaglio et aL 2005b. Nevertheless we note that 



this later discrepancy with lErb et alJ (l2006h results may be un- 
derstood: according to the downsizing scenario, the evolution of 
the most massive galaxies plotted here should be actually smaller 
betwee n z = 2 and z = 1 than betwe en z = 1 and z = 0. Among 
others, iPerez-Gonzalez et al.l (|2008|) have quantified that galax- 
ies below log(M*/M©) = 11 .5 have formed half of th eir stars at 
z < 1. The main reason of a possible overestimate of lErb et"an 
(120061) metallicities is probably statistical variation effects. They 
have indeed based their results on stacked spectra of very few 
galaxies. The effect of the selection function is therefore diffi- 
cult to analyze. 

The type of galaxies observed by lErb et~an (l200d) at z ~ 2, 
which are active galaxies, may also have later evolved to "red 
and dead" passive galaxies, and got hi gher metalliciti e s than 
the ones actually observed by us or bv iTremonti et al.l (|2004|) 
at lower redshifts. Such dead galaxies would unfortunately not 
satisfy any more the selection function of any work based on 
emission-line measurements, and would not be observed. 



M(B) 




M(B) 



Figure 16. This plot shows, for galaxies with similar stellar 
masses, their evolution when redshift increases: for metallicity 
in the mass-metallicity plane (top-left), for metallicity in the 
luminosity-metallicity plane (top-right), and for luminosity in 
the mass-luminosity plane (bottom-right). 

Table 9. Evolution of the mass-to-light ratio for the wide and 
deep samples. This evolution is computed as an absolute rest- 
frame B-band magnitude evolution at constant stellar mass, from 
the comparison of the luminosity-metallicity (see Fig. [12] and 
Table|6| and mass-metalUcity (see Fig.[T3]and Table|7]i relations. 



redshift range 



AMab(S) 
wide deep 



0.0<z<0.5 
0.5 <z<0.7 
0.7 < z < 0.9 



-0.06 0.65 
-0.61 -0.10 
-1.19 -0.45 



6.4. Evolution of the mass-to-light ratio 

The metallicity evolution on the wide sample is very similar 
to what has been found by iLamareille et al.l ( l2006al) on the 
luminosity-metalUcity relation, after having applied a correction 



for luminosity evolution jllbert et alj|2005l) . This tells us that the 
stronger evolution of the luminosity-metallicity relation, com- 
pared to the mass-metallicity relation, is effectively due to an 
additional luminosity evolution. This luminosity evolution has 
to be understood for galaxies with similar masses, which means 



we can measure the evolution of the mass-to-light ratio in galax- 
ies using metallicity as a pivot. 

Fig. [16] shows schematically the global evolution of metal- 
licity in the luminosity-metallicity plane (top-right). When there 
is both an evolution of the metallicity in the mass-metallicity 
plane (top-left), and an evolution of the luminosity in the mass- 
luminosity plane (bottom-right), galaxies being more luminous 
at a given mass, the result is a stronger evolution of the metallic- 
ity in the luminosity-metallicity plane. 

The luminosity evolution at constant mass of our galaxies 
can be calculated using the following formula: 

AMab(B) =a-'x (Alog(0/H)'*^ - Alog(0/H)^) (6) 

where a is the slope of the luminosity-metallicity relation, i.e. 
—0.31 dex/mag. The results are given in Table |9] The positive 
evolution for the lowest redshift bin in the deep sample is due 
to the incompleteness of the N2 calibration: at a given metallic- 
ity, our sample is biased towards higher luminosity. In the other 
redshift ranges and the whole wide sample, we see that galax- 
ies at higher redshifts have lower mass-to-light ratios than today, 
and that this evolution is more significant for massive galaxies 
(i.e. the wide sample compared to the deep). These two results 
are in good agreem ent with the general scenario of downsizing 
dCowie et al .119961) for the evolution of the star formation rates in 
galaxies. Less massive galaxies show a smaller evolution in their 
star formation rate at z < 1 since they are still actively forming 
stars. 

It will be further analyzed in a subsequent paper of this se- 
ries. 

6.5. Derived star formation rates 

We now discuss the evolution of metallicity at constant stellar 
mass, in terms of star formation rates. 

Assuming constant star formation rates (SFR) and the 
closed-box model, one can calculate a stellar mass evolution and 
relates it to a metallicity evolution using the following equations: 

M^{t)=M^{Q)+SFRxt (7) 
M*(f)+M,,(f)-Mtot (8) 

where M^,, Mg and Mtot are respectively the stellar mass, the gas 
mass, and the total baryonic mass of the galaxy (which remain 
constant in the closed-box model); and where Z and y are respec- 
tively the metallicity and the true yield (which depends only on 
the stellar initial mass function). 

We know the m ass-metallicity relation at redshift z = 
(iTremonti et alJl2004l) . and the Table|7]gives the metallicity evo- 
lution for various stellar masses and redshifts. Thus, we can re- 
vert Eq. |7] [8] and |9] to derive the star formation rate which ex- 
plains the observed values. The results are quoted in Table [TO] 
(case fl). They are cal culated for a given c osmology, and for a 
true yield y = 0.0104 (ITremonti et alJl2004 . 

A better modeling can be performed by taking into account 
the time evolution of the star formation rates in galaxies. We 
may for example assume an exponentially decreasing star for- 
mation rate, i.e. SFR{t) ~ SFR{Q) x exp(— f/r), where T gives 
the characteristic e-folding time of the galaxy. The exponentially 
decreasing law is a good choice when analyzing a population of 



galaxies, with respect to the global cosmic evolution of star for- 
mation rate. Eq. |2]is then replaced by the following formula: 

M^{t)^M^{0)+SFR{0)xTx (^l-g-'/^) (10) 

Table [TOl (case b) gives the r esults obtained when assuming 
the relation given by Eq. 12 of ISavaglio et al 1 (I2005h between 
the e-folding time and the total baryonic mass of the galaxy. 
Comparing cases a and b, we see that the results are not dra- 
matically affected by the assumption of a decreasing SFR. The 
star formations activities in case b are systematically higher than 
in case a, which is expected as with a decreasing SFR a higher 
initial value is needed to explain the same metallicity evolution, 
as compared to a constant SFR. 

We note also that the derived star formation activities in the 
wide and deep samples are very close, despite their different 
ranges in stellar masses. In case a, star formation activities in 
the deep sample are nevertheless higher than in the wide sam- 
ple, which is expected as the wide sample spans higher stellar 
masses. Conversely, in case b, we find higher star formation 
activities for more massive objects. Th is comes from the as- 
sumption, made by ISavaglio et al ] (l2005h and used in our equa- 
tions, that less massive galaxies have a longer e-folding times. 
Consequently the less massive galaxies of the deep sample, 
which are assigned long e-folding times ans thus can be approx- 
imated as constant SFR, do not vary much from case a to case b 
while the more massive galaxies of the wide sample, which are 
assigned short e-folding times, are more affected by the decreas- 
ing SFR hypothesis. 

Preliminary results on the evolution of star formation activi- 
ties of VVDS g alaxies, which w ill be presented in a subsequent 
paper (see also IWalcher et alJ llOOSj), give at redshift z ~ 0.6 
(our best sampled redshift bin) the following numbers: « 10^^^ 
yr^^ for galaxies at « 10^'^ solar masses (mean mass of the deep 
sample), and « lO^"*'^ yr^' for galaxies at w lO'" solar masses 
(mean mass the wide sample). 

Assuming the closed-box hypothesis, our results in the wide 
sample are in fairly good agreement with the observed star for- 
mation activities at similar masses, this agreement being better 
in case b. But in the deep sample, the inferred star formation ac- 
tivities which explain our observed metallicity evolution seem 
significantly underestimated by a factor « 5. Although the vary- 
ing e-folding times hypothesis is in agreement with the slower 
evolution of the mass-to-light ratio found in the deep sample in 
previous section, it fails in providing a good agreement between 
our inferred star formation activities in the deep sample and ob- 
served values at similar masses. Moreover this hypothesis make 
more massive galaxies forming stars more actively, which is in 
contradiction with observations. 

We therefore conclude that the closed-box hypothesis is not 
valid to explain our observed metallicity evolutions (assuming 
that the preliminary results mentioned above are confirmed). 

7. Conclusion 

We have calculated rest-frame luminosities, stellar masses 
and gas-phase oxygen abundances of a statistically significant 
sample of star-forming galaxies, selected from the VIMOS 
VLT Deep Survey. This has allowed us to derive luminosity- 
metallicity and mass-metallicity relations in various redshift 
ranges up to z ~ 0.9. These relations have also been derived in 
two sub-samples: the deep sample which extends to lower ob- 
served luminosities thanks to a deeper magnitude selection, and 
the wide sample which extends to higher luminosities thanks to 



Table 10. Derived star formation activities needed to explain the observed values quoted in Table |71 assuming the closed-box 
model, and h = 0.7, D.^ = 0.3 and Q.a = 0.7 cosmology. The values are given in the wide and deep samples, and for two different 
assumptions: a) constant star formation rate, b) exponentially declining star formation rate (see text). In case b, each value is the 
instantaneous star formation rate at the specified redshift (see Eq.fTOli. 



mean redshift 




wide 




log(SFR/M*) 
deep 




a 


b 


a 


b 


0.3 


-11.06±0.7 


-10.93±0.7 


-11.05±0.9 


-10.99±0.9 


0.6 


-10.64 ±0.5 


-10.36±0.5 


-10.51 ±0.4 


-10.40 ±0.5 


0.77 


-10.73±0.4 


-10.27±0.5 


-10.65±0.5 


-10.53±0.5 



a larger solid angle. The selection function have been taken into 
account in order to define volume- and mass-limited samples. 

For both the wide and the deep samples, we have studied 
as a function of redshift: the evolution in metallicity at constant 
luminosity, the evolution in metallicity at constant stellar mass, 
the associated evolution in luminosity at constant stellar mass, 
and the evolution of the slope of the mass-metallicity relation. 

Additionally, we have also: 

- Measured the fraction of star-forming galaxies and AGNs as 
a function of redshift. 

- Calculated the difference in terms of stellar masses between 
the old BC03 and new CB07 stellar population models. 

- Studied the effect on the computation of stellar mass of the 
use of spectral indices rather than only photometric points. 

- Checked that gas-phase oxygen abundance measurements in 
low and high redshift ranges are comparable even if calcu- 
lated with different methods. 

- Studied the marginal effect of candidate star-forming 
galaxies and candidate AGNs on the derived luminosity- 
metalhcity or mass-metallicity relations. 

In the wide sample, the mean evolution of metallicity at con- 
stant luminosity, and of luminosity at constant stellar mass, are 
in good agreement w i th prev i ous studies by e . g.,Lamareille et al.l 
(l2006ah : lllbert et all (l2005h : iHammer et all OOOSh among oth- 
ers. 

By doing the assumption that the slope of the luminosity- 
metallicity relation, or the shape of the mass-metallicity rela- 
tion, remain constant with redshift, we finally found different 
metallicity evolutions for the wide and deep samples. The wide 
and deep samples span different ranges in masses, which shows 
that the assumption is wrong and that this slope or shape actu- 
ally evolves with redshift. We have found that the most massive 
galaxies show the strongest metallicity evolution. At z ~ 0.77, 
galaxies at 10^^ solar masses have —0.18 dex lower metallici- 
ties than galaxies of similar masses in the local universe, while 
galaxies at lO"^'^ solar masses have —0.28 dex lower metallici- 
ties. 

We have then studied the mass-metallicity relation on co- 
added data points by bins of stellar masses. The inferred slopes 
show an evolution to a flatter mass-metallicity relation at z ~ 1 as 
compared to the local universe, which confirms that higher mass 
galaxies show a stronger metallicity evolution during this period. 
The slope evolves by —0.25 dex/decade by unit of redshift in the 
range < z < 0.7. 

Moreover, we have inferred from our observations the star 
formation activities which explain this evolution in metallicity. 
To do so, we have assumed the closed-box hypothesis. We have 
found that the observed metallicity evolution would be explained 
by similar star formation activities at any mass, which is in con- 
tradiction with the well-known correlation between stellar mass 



and star formation activity. Assuming preliminary observations 
of the star formation activities in VVDS galaxies, we conclude 
that our results are underestimated by a factor « 5 in the deep 
sample, while their are in fair agreement in the wide sample. 

These results make a strong evidence against the closed- 
box model. Indeed, we know from previous studies that the star 
formation activity of galaxies decreases while stellar mass in- 
creases. The only way to explain the smaller metallicity evolu- 
tion of the less massive galaxies is thus to replace in Eq. |9] the 
true yield by a smaller effective yield. Assuming a smaller effec- 
tive yield allows galaxies to show smaller metallicity evolution, 
even with a stronger star formation activity. 

The smaller effective yield of less massive galaxies can be 
understood in the "open-closed" model. In this model, galax- 
ies with small stellar masses evolve like open-boxes: most of 
the metals produced during star formation are ejected in the in- 
tergalactic medium by stellar winds and supernovae feedback, 
the effective yield is thus very small. Conversely, galaxies with 
high stellar masses evolve in the open-closed model like closed- 
boxes: the metals are retained in the galaxy thanks to a high 
gravitational potential, the effective yield is thus close to the true 
yield. 

The dependence of the effective yield with the gravitational 
potential of the galaxies, itself related to thei r total baryonic 
mass, has been already shown by lTremonti et al.l (l2004l) . It is nat- 
urally explained in the hierarchical galaxy formation scenario: in 
contrary to the closed-box model, the total baryonic mass does 
not remain constant and increases with time, together with the 
stellar mass, thanks to galaxy merging and accretion. This rela- 
tion between the stellar mass and the total baryonic mass is at the 
origin of the mass-metallicity relation in the open-closed model, 
and is clearly supported by our data. 

In the open-closed model, the smaller metallicity evolution 
of less massive is naturally explained despite their stronger star 
formation activities. 

We emphasize that Eq. |7] [8] and [TO] are not valid any more 
in the open-closed model. It is not enough to replace, in Eq. |9] 
the true yield by the effective yield in order to infer the right 
star formation rates which explain the observed metallicity evo- 
lution. These equations must be modified to take also the dif- 
ference between the true and the effective yield into account. 
The way to modify these equations depends on the model used 
to explain smaller effective yields (e.g. gas loss). Replacement 



equations can be found in other studies (e.g.|Erb et al 
■ 12 



eplacement 
i1(2006tl&bl 



l2008l:lDalcantonl2007HFinlator & Davel2008ir We also refer the 
reader to subsequent papers of this series, for a more detailed 
discussion of the relation between the present results and the 
evolution of the star formati on rates of galaxies, a nd for compar- 
isons with simulations (e.g. lDe Lucia et a 02004 . 
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